On the use of sensitivity tests in seismic tomography

ACKNOWLEDGEMENTS This work was partly supported by ARC Discovery Project DP120103673 and by the Research Council of Norway through its Centres of Excellence funding scheme, project number 223272. We thank Maximilliano Bezada and an anonymous referee for constructive comments which improved the original version of the manuscript. We also thank the Editor, A. Morelli, for providing additional helpful comments.Peer reviewedPublisher PD

Global correlations between maximum magnitudes of subduction zone interface thrust earthquakes and physical parameters of subduction zones

Peer reviewedPublisher PD

Ambient noise tomography reveals basalt and sub-basalt velocity structure beneath the Faroe Islands, North Atlantic

© 2017 Elsevier B.V. Ambient noise tomography is applied to seismic data recorded by a portable array of seismographs deployed throughout the Faroe Islands in an effort to illuminate basalt sequences of the North Atlantic Igneous Province, as well as underlying sedimentary layers and Precambrian basement. Rayleigh wave empirical Green's functions between all station pairs are extracted from the data via cross-correlation of long-term recordings, with phase weighted stacking implemented to boost signal-to-noise ratio. Dispersion analysis is applied to extract inter-station group travel-times in the period range 0.5–15 s, followed by inversion for period-dependent group velocity maps. Subsequent inversion for 3-D shear wave velocity reveals the presence of significant lateral heterogeneity (up to 25%) in the crust. Main features of the final model include: (i) a near-surface low velocity layer, interpreted to be the Malinstindur Formation, which comprises subaerial compound lava flows with a weathered upper surface; (ii) a sharp velocity increase at the base of the Malinstindur Formation, which may mark a transition to the underlying Beinisvørð Formation, a thick laterally extensive layer of subaerial basalt sheet lobes; (iii) a low velocity layer at 2.5–7.0 km depth beneath the Beinisvørð Formation, which is consistent with hyaloclastites of the Lopra Formation; (iv) an upper basement layer between depths of 5–9 km and characterized by S wave velocities of approximately 3.2 km/s, consistent with low-grade metamorphosed sedimentary rocks; (v) a high velocity basement, with S wave velocities in excess of 3.6 km/s. This likely reflects the presence of a crystalline mid-lower crust of Archaean continental origin. Compared to previous interpretations of the geological structure beneath the Faroe Islands, our new results point to a more structurally complex and laterally heterogeneous crust, and provide constraints which may help to understand how continental fragments are rifted from the margins of newly forming ocean basins

Seismic imaging at the cross-roads: Active, passive, exploration and solid Earth

Science has grown from our need to understand the world around us. Seismology is no different, with earthquakes and their destructive effect on society providing the motivation to understand the Earth's seismic wavefield. The question of when seismology as a science really began is an interesting one, but it is unlikely that there will ever be a universally agreed-upon date, partly because of the incompleteness of the historical record, and partly because the definition of what constitutes science varies from person to person. For instance, one could regard 1889 as the true birth of seismology, because that is when the first distant earthquake was detected by an instrument; in this case Ernst von Rebeur-Paschwitz detected an earthquake in Japan using a pendulum in Potsdam, Germany (Ben-Menahem, 1995). However, even the birth of instrumental seismology could be contested; the so-called Zhang Heng directional “seismoscope” (detects ground motion but not as a function of time) was invented in 132 CE (Rui and Yan-xiang, 2006), and is said to have detected a four-hundred mile distant earthquake which was not felt at the location of the instrument Needham, 1959, Dewey and Byerly, 1969. Prior to instrumental seismology, observations of earthquakes were not uncommon; for instance, Aristotle provided a classification of earthquakes based on the nature of observed ground motion (Ben-Menahem, 1995)

The mechanisms underpinning Cenozoic intraplate volcanism in eastern Australia: Insights from seismic tomography and geodynamic modeling

Cenozoic intraplate volcanism is widespread throughout much of eastern Australia and manifests as both age-progressive volcanic tracks and non-age-progressive lava fields. Various mechanisms have been invoked to explain the origin and distribution of the volcanism, but a broad consensus remains elusive. We use results from seismic tomography to demonstrate a clear link between lithospheric thickness and the occurrence, composition, and volume of volcanic outcrop. Furthermore, we find that non-age-progressive lava fields overlie significant cavities in the base of the lithosphere. Based on numerical simulations of mantle flow, we show that these cavities generate vigorous mantle upwellings, which likely promote decompression melting. However, due to the intermittent nature of the lava field volcanics over the last 50 Ma, it is probable that transient mechanisms also operate to induce or enhance melting. In the case of the Newer Volcanics Province, the passage of a nearby plume appears to be a likely candidate. Our results demonstrate why detailed 3-D variations in lithospheric thickness, plate motion, and transient sources of mantle heterogeneity need to be considered when studying the origin of non age-progressive volcanism in continental interiors

Evidence of micro-continent entrainment during crustal accretion

Peer reviewedPublisher PD

Inherited crustal deformation along the East Gondwana margin revealed by seismic anisotropy tomography

Acknowledgments We thank Mallory Young for providing phase velocity measurements in mainland Australia and Tasmania. Robert Musgrave is thanked for making available his tilt-filtered magnetic intensity map. In the short term, data may be made available by contacting the authors (S.P. or N.R.). A new database of passive seismic data recorded in Australia is planned as part of a national geophysics data facility for easy access download. Details on the status of this database may be obtained from the authors (S.P., N.R., or A.M.R.). There are no restrictions on access for noncommercial use. Commercial users should seek written permission from the authors (S.P. or N.R.). Ross Cayley publishes with the permission of the Director of the Geological Survey of Victoria.Peer reviewedPublisher PD

Insights into the structure and dynamics of the upper mantle beneath Bass Strait, southeast Australia, using shear wave splitting

© 2019 Elsevier B.V. We investigate the structure of the upper mantle using teleseismic shear wave splitting measurements obtained at 32 broadband seismic stations located in Bass Strait and the surrounding region of southeast Australia. Our dataset includes ∼366 individual splitting measurements from SKS and SKKS phases. The pattern of seismic anisotropy from shear wave splitting analysis beneath the study area is complex and does not always correlate with magnetic lineaments or current N-S absolute plate motion. In the eastern Lachlan Fold Belt, fast shear waves are polarized parallel to the structural trend (∼N25E). Further south, fast shear wave polarization directions trend on average N25–75E from the Western Tasmania Terrane through Bass Strait to southern Victoria, which is consistent with the presence of an exotic Precambrian microcontinent in this region as previously postulated. Stations located on and around the Neogene-Quaternary Newer Volcanics Province in southern Victoria display sizeable delay times (∼2.7 s). These values are among the largest in the world and hence require either an unusually large intrinsic anisotropy frozen within the lithosphere, or a contribution from both the lithospheric and asthenospheric mantle. In the Eastern Tasmania Terrane, nearly all observed fast directions are approximately NW-SE. Although part of our data set strongly favours anisotropy originating from “fabric” frozen in the lithospheric mantle, a contribution from the asthenospheric flow related to the present day plate motion is also required to explain the observed splitting parameters. We suggest that deviation of asthenospheric mantle flow around lithospheric roots could be occurring, and so variations in anisotropy related to mantle flow may be expected. Alternatively, the pattern of fast polarisation orientations observed around Bass Strait may be consistent with radial mantle flow associated with a plume linked to the recently discovered Cosgrove volcanic track. However, it is difficult to characterise the relative contributions to the observed splitting from the lithospheric vs. asthenospheric upper mantle due to poor backazimuthal coverage of the data

Automatic differentiation in geophysical inverse problems

Automatic differentiation (AD) is the technique whereby output variables of a computer code evaluating any complicated function (e.g. the solution to a differential equation) can be differentiated with respect to the input variables. Often AD tools take the form of source to source translators and produce computer code without the need for deriving and hand coding of explicit mathematical formulae by the user. The power of AD lies in the fact that it combines the generality of finite difference techniques and the accuracy and efficiency of analytical derivatives, while at the same time eliminating 'human' coding errors. It also provides the possibility of accurate, efficient derivative calculation from complex 'forward' codes where no analytical derivatives are possible and finite difference techniques are too cumbersome. AD is already having a major impact in areas such as optimization, meteorology and oceanography. Similarly it has considerable potential for use in non-linear inverse problems in geophysics where linearization is desirable, or for sensitivity analysis of large numerical simulation codes, for example, wave propagation and geodynamic modelling. At present, however, AD tools appear to be little used in the geosciences. Here we report on experiments using a state of the art AD tool to perform source to source code translation in a range of geoscience problems. These include calculating derivatives for Gibbs free energy minimization, seismic receiver function inversion, and seismic ray tracing. Issues of accuracy and efficiency are discussed

P and S wave travel time tomography of the SE Asia-Australia collision zone

© 2019 Elsevier B.V. The southeast (SE)Asia - Australia collision zone is one of the most tectonically active and seismogenic regions in the world. Here, we present new 3-D P- and S-wave velocity models of the crust and upper mantle by applying regional earthquake travel-time tomography to global catalogue data. We first re-locate earthquakes provided by the standard ISC-Reviewed and ISC-EHB catalogues using a non-linear oct-tree scheme. A machine learning algorithm that clusters earthquakes depending on their spatiotemporal density was then applied to significantly improve the consistency of travel-time picks. We used the Fast Marching Tomography software package to retrieve 3-D velocity and interface structures from starting 1-D velocity and Moho models. Synthetic resolution and sensitivity tests demonstrate that the final models are robust, with P-wave speed variations (~130 km horizontal resolution)generally recovered more robustly than S-wave speed variations (~220 km horizontal resolution). The retrieved crust and mantle anomalies offer a new perspective on the broad-scale tectonic setting and underlying mantle architecture of SE Asia. While we observe clear evidence of subducted slabs as high velocity anomalies penetrating into the mantle along the Sunda arc, Banda arc and Halmahera arc, we also see evidence for slab gaps or holes in the vicinity of east Java. In the Banda arc, we image the slab as a single curved subduction zone. Furthermore, a high-velocity region in the mantle lithosphere connects northern Australia with Timor and West Papua. The S-wave model shows broad-scale features similar to those of the P-wave model, with mantle earthquakes generally distributed within high-velocity slabs. The high velocity mantle connection between northern Australia and the eastern margin of the Sunda arc is also present in the S-wave model. While the S-wave model has a lower resolution than the P-wave model due to the availability of fewer paths, it nonetheless provides new and complementary insights into the structure of the upper mantle beneath southeast Asia